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Abstract 

A method is developed to estimate the parameters of a Levy copula of a discretely 
observed bivariate compound Poisson process without knowledge of common shocks. 
The method is tested in a small sample simulation study. Also, the method is applied 
to a real data set and a goodness of fit test is developed. With the methodology of 
this work, the Levy copula becomes a realistic tool of the advanced measurement 
approach of operational risk. 

o 
q 

^ 1 Introduction 

The multivariate compound Poisson process is an intuitively appealing and natural model 
for operational risk and insurance claim modelling. The model is intuitively appealing 
because dependencies between different loss categories are caused by common shocks 
that apply to multiple loss categories simultaneously. For example, in operational risk 
modelling, failure of an IT system is a common shock that causes losses in multiple 
lines of business. The multivariate compound Poisson process is a natural model for the 
following two reasons. First, as a Levy process, it is easily applied to any time horizon 
of interest. Second, because a redesign of loss categories results in a loss process that is 
again multivariate compound Poisson, the nature of the model does not depend on the 



i 



level of granularity (Bocker and Kliippelberg, 2008). 
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A multivariate compound Poisson process can be specified in terms of univariate com- 
pound Poisson processes and a Levy copula (Cont and Tankov, 2004). In essence, a Levy 
copula provides the relationship between the Levy measure of a multivariate Levy process 
and the Levy measures of its marginal processes. The Levy copula allows for a parsimo- 
nious bottom-up modelling with compound Poisson processes. In case of two loss cate- 
gories, for example, parameterization with a Clayton Levy copula requires two marginal 
frequencies, two marginal jump size distributions and one Levy copula parameter. In con- 
trast, parameterization without a Levy copula requires three frequencies (corresponding 
to losses of the first category only, losses of the second category only and common shocks 
that apply to both categories), two univariate jump size distributions (corresponding to 
losses of one of the two categories only) and, finally, one bivariate jump size distribution 
(corresponding to the common shocks). 

The parameters of a Levy copula of a multivariate compound Poisson process can 
be estimated if the process is either observed continuously (such that common shocks 
can be identified) or observed discretely with knowledge about all jump sizes and the 
common shocks. These two cases have been studied by Esmaeili and Kliippelberg (2010) 
for a bivariate compound Poisson process (the continuous observation is mimicked in a 
simulation study, while the discrete observation corresponds to a real data set of insurance 
claims). The objective of this work is to develop a method to estimate the parameters 
of a Levy copula of a bivariate compound Poisson process in case the process is observed 
discretely with knowledge about all jump sizes, but without knowledge of which jumps 
stem from common shocks. This situation is relevant to operational risk modelling in 
which all material losses are registered, but common shocks are typically unknown. With 
the methodology developed here, the Levy copula becomes a realistic tool of the advanced 
measurement approach of operational risk. 

The outline of this paper is as follows. In Section |2j we discuss the bivariate compound 
Poisson process in terms of its common shock representation and Levy measure. This 
prepares the ground for the two-dimensional Levy copula of Section [3] In Section |4| the 
new estimation method of the discretely observed bivariate compound Poisson process is 
presented. The method is tested in a simulation study in Section [5] and applied to a real 
data set in Section [6j In Section |6j we also develop a goodness of fit test for the Levy 
copula. Finally, we conclude in Section [7j 
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2 The bivariate compound Poisson process 



A bivariate compound Poisson process S = (S(t)) t >o = ((Si(t), S 2 (t))) t >o is defined on a 
filtered probability space (ft, J 7 , P) as 

N(t) 

s(*) = X>, C 1 ) 

1=1 

where N = (N(t)) t >o is a Poisson process with intensity A > and Y = (Yi) i( z N + is a 
sequence of iid 2-dimensional random vectors. The process N and the sequence Y are 
statistically independent. By construction, given any t > s, the increment S(t) — S(s) is 
independent of T(s) and has the same distribution as S(t—s). The probability distribution 
of Yi is such that F(Yi = 0) = 0, which means that a jump of N almost surely manifests 
itself in a jump of at least one of the components of S. 

2.1 The Levy-Ito decomposition 

The Levy-Ito decomposition of S(t) takes the form 

S(t) =11 xJ{ds,dx), (2) 



where J is the Poisson random measure. With the help of the Levy-Ito decomposition, 
we find that S(t) has common shock representation 



where 



and 



S 1 (t) = St(t) + S l }(t) and S 2 (t) = S 2 L (t) + S l 2 \t), (3) 



Sf(t)= I [ xiJ(rfs,rfxi x {0}), (4) 
Jo Jr\{o} 

S 2 L {t)= [ [ x 2 J{ds, {0} x dx 2 ) (5) 
Jo Jr\{o} 

S ]l (t)=.(sl(t),S 2 l (t))= [ f [x u x 2 )J{ds,dxi x dx 2 ). (6) 

Jo J(R\{0}) 2 



The processes = (S , j L (t)) t > , = (S 2 (t)) t > and S" = (S^(t)) t > do not jump 



simultaneously and are statistically independent (Sato, 1999). The processes Sf and 5*2 
are called the independent parts of S. Conversely, the process S"" is called the dependent 
part of S and corresponds to the common shocks. 
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2.2 The Levy-Khinchin representation and the Levy measure 

The Levy-Khinchin representation of the characteristic function of S(t) can be determined 
from Eq. ^ with the exponential formula for Poisson random measures ( Cont and Tankov 



2004). The representation takes the form 



E [ e iu - s ®] = exp t I (e iu - x - l) u(dx) 



(7) 



where u G IR 2 and the Levy measure 

u(A) = XP(Yi G A) 



E[J((0,t],A)} 



t 



gives the expected number of jumps per unit of time in each Borel set A of M 2 . The 
processes Si and S2 are independent if and only if the support of v is contained in the 
set {x G 



52 . 



: X1X2 = 0} (Cont and Tankov, 2004). In this case, Si and S 2 do not jump 
simultaneously almost surely and the Levy-Khinchin representation factorizes as 



E [. 



,iu 1 Si(t)+iu 2 S 2 {t) 



] = E [ e iuiSl W] E [ e in252{t) ] 



(9) 



On the other hand, the processes Si and S 2 are defined to be comonotonic if their jump 



sizes ASi and AS 2 , respectively, are elements of an increasing set S of R , see (Cont and 



Tankov, 2004). Any two elements (xi,x 2 ) and (1/1,2/2) of S satisfy Xi > yi or X{ < yi for 
all i = 1,2. An example of an increasing set is {x G IR 2 : Xi = x 2 }. The requirement 
(ASi, AS 2 ) G S means that by observing one of the processes Si or S 2 , the other process 
can be constructed exactly with a positive dependence. In case of comonotonic S% and 
S2, the Levy measure is concentrated on S. 

In terms of Si(t), S 2 (t) and S 1 ", Eq. (TH) takes the form 



E [ 



z ,iu 1 S 1 {t)+iu 2 S 2 {t) 



] =E 



E 



E 



(10) 



where we have used that S^, and S" are independent. The Levy-Khinchin represen- 
tation of the characteristic functions of S±(t), S 2 (t) and S"(t) can be determined from 
their Levy-Ito decompositions in the same way as Eq. ^ is determined from Eq. ^ . The 
Levy measures of S± and S^ are given by, respectively, 



v±(B) = u(B x {0}) and u£(B) = z/({0} x B), 
where B is a Borel set of R. The Levy measure of S" takes the form 



J (A) = u(A) - v{Ai x {0}) - v({0} x A 2 ), 



(12) 
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where the sets A\ and A 2 are defined as 

At = {x G M 2 : (zi, 0) G A} and A 2 = {i 6 I 2 : (0, x 2 ) G A}. (13) 



To conclude our discussion of the bivariate compound Poisson process S, we consider its 
components for a Levy measure v that is not necessarily concentrated on {i G K 2 : x x x 2 = 



0} or an increasing set S. The process Si is compound Poisson (Cont and Tankov, 2004) 



and by setting u = («i,0) in Eq. (10), we find that the characteristic function of S\(t) 
takes the form 



E[e iUlSl{t) ] =exp t [ (e iUlXl -l)vj-(dxi) exp t [ (e™ - l) J {dx x x dx 2 ) 

i Jr ' J L Jr 2 



exp 



t / (e tUlXl - 1) {vi{dxi) + J{dxi x (- 



-oo, oo 



))} 



(14) 



From Eq. ( 14 ) , it follows that the Levy measure of Si is given by 



(15) 



If the measure v is concentrated on {x G 1R 2 : XiXi = 0}, then Vi(B) = fi(B) and if it is 
concentrated on (IR\{0}) 2 (such as on an increasing set S), then v x {B) = J(B xl). In 



general, ui is a combination of z/f and cf. Eq. (15). In the same way, 5*2 is compound 
Poisson with Levy measure 



v 2 {B) = v^(B) + i/H(R x B). 



(16) 



3 The Levy copula 

We consider a bivariate compound Poisson process with positive jumps. This means that 
the Levy measure is concentrated on [0, oo) 2 \{0} rather than on M 2 \{0}. The assumption 
of positive jumps is reasonable in the context of operational risk modelling and restricts 
our discussion of Levy copulas to positive Levy copulas. 



3.1 The positive Levy copula and Sklar's theorem 

A two-dimensional positive Levy copula C : [0, oo] 2 — > [0, oo] is a 2-increasing grounded 
function with uniform margins. The 2-increasing property means that for any (^1,^2) G 
[0, oo] 2 and (f 1,1*2) G [0, oo] 2 with Ui > v x and u 2 > v 2 , we have 

C(ui,u 2 ) - C(ui,v 2 ) -C(vi,u 2 ) + C(vi,v 2 ) > 0. (17) 
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The grounded property means that C{u\,u 2 ) — if u\ = and/or u 2 = 0. Finally, 
margins are defined as C(ui, oo) and C(oo,m 2 ), and the positive Levy copula is such that 
C(ui, oo) = ui and C(oo, u 2 ) = u 2 . 

In the same way as a distributional copula connects marginal distribution functions 
to a joint distribution function, the Levy copula connects marginal tail integrals to a tail 
integral. For a two-dimensional Levy process with positive jumps, the tail integral U is 
defined as 

{0 if x\ = oo and/or x 2 = oo 

u([ Xl ,oo) x [x 2 ,oo)) if (x\,x 2 ) G [0,oo) 2 \{0} (18) 
oo if (xi, x 2 ) = 

and the marginal tail integrals are defined as 

Ui(x x ) = U{x x ,U) and U 2 (x 2 ) = U(0, x 2 ). (19) 



The following theorem due to Cont and Tankov (2004) is a version of Sklar's theorem for 
Levy copulas. 

Theorem 1 Let (S\,S 2 ) be a two-dimensional Levy process with positive jumps having 
tail integral U and marginal tail integrals U\ and U 2 . There exists a two-dimensional 
positive Levy copula C which characterizes the dependence structure of (Si,S 2 ), that is, 
for all X\, x 2 G [0, oo], 

U(x 1 ,x 2 )=C(U 1 (x 1 ),U 2 {x 2 )). (20) 

IfU\ and U 2 are continuous, this Levy copula is unique. Otherwise it is unique on Ranf/i x 
Ran L^. 

Conversely, let Si and S 2 be two one- dimensional Levy processes with positive jumps 
having tail integrals U\ and U 2 and let C be a two-dimensional positive Levy copula. Then 
there exsists a two-dimensional Levy process with Levy copula C and marginal tail integrals 



U\ and U 2 . Its tail integral is given by Eq. (20) 



The definition of the tail integral and its marginal tail integrals imply that C/i(0) = oo 
and f/ 2 (0) = oo. The singularity at zero is necessary to correctly account for jumps of the 
independent parts and S^- Consider, for example, on the one hand 

C(£/i(xi),£/ 2 (0)) = C(£/i(xi),oo) = C/i(xi) = u([x u oo)x [0,oo)) for < xi < oo (21) 

and, on the other hand 



lim C(L r i(xi), U 2 (x 2 )) = lim U(x\,x 2 ) = lim u([xi, oo) x [x 2 , oo)) for < x\ < oo. 

X2^0 X2-10 X2^0 



(22) 
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The difference 

v([xi, oo) x [0, oo)) — lim u([xi, oo) x [x 2 , oo)) = v([xi, oo) x {0}) (23) 

£24-0 



for < X\ < oo between the tail integrals of Eqs. (21) and (22) corresponds to the tail 
integral of S^. If U 2 (0) would have been defined as lim X2 ^ ^2(^2), the difference of the 
tail integrals vanishes and Sf~ does not jump almost surely. 

3.2 Construction of positive Levy copulas 

In case of independent Si and S 2 , the Levy measure is concentrated on the set {x G 
[0, oo) 2 : x\x 2 = 0} and the tail integral takes the form 

U (xi, x 2 ) = £/i (£1)1^=0 + U 2 {x 2 )l xi=0 . (24) 



With the help of Eq. (20), we find that the independence Levy copula C± is given by 

C±(ui,u 2 ) = Mil U2=00 + u 2 l (25) 

In case of comonotonic S\ and S 2 , the Levy measure is concentrated on an increasing set 
S and the tail integral takes the form 

U{x u x 2 ) =mm{U 1 {x l ),U 2 {x 2 )), (26) 

which implies that the comonotonic Levy copula Cu is given by 

C\\ (ui, u 2 ) = min(ui, u 2 ). (27) 

A Levy copula C with a dependence that is between C± and Cn can be constructed in 
several ways, such as by an approach similar to the construction of Archimedean distri- 



butional copulas (Cont and Tankov, 2004). Given a strictly decreasing convex function 
(f> : [0, 00] — > [0, 00] such that 0(0) = 00 and <fi(oo) = 0, a positive two-dimensional 
Archimedean Levy copula is defined as 

C( Ul ,u 2 ) =0- 1 (0(m 1 ) + 0(m 2 )). (28) 

For 4>{u) = u~ 5 with 5 > 0, one obtains the Clayton Levy copula 

C(ui,u 2 ) = (uf + u 2 5 )- 1/5 ■ (29) 

The Clayton Levy copula includes the independence Levy copula C± and the comonotonic 
Levy copula C\\ in the limits 5 J, and 5 — > 00, respectively. 
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3.3 Dependence implied by the Levy copula 

The bivariate compound Poisson process is fully determined by the Levy measures , 
i>2 and v». Given a Levy copula C, these measures can be expressed in terms of the Levy 
measures ui and z/ 2 . All Levy measures are defined in the same way as in Section [2j but 
now v is concentrated on [0, oo) 2 \{0}. The relation between the measures , z/^ and 
on the one hand, and the measures v x and z/ 2 on the other, is established in Appendix 
|A~} In this Section, for later purposes, the relation between the measures is expressed in 
terms of frequencies and jump size distributions. 

The frequencies A^ and A 2 of, respectively, S± and S 1 ^ are given by 

= ^((0,00)) = Ai - A l! and A^ = ^((0, 00)) = A 2 - A 11 , (30) 

where 

Ai = ^((0,oo)), A 2 = z/ 2 ((0,oo)) and A 11 = J ((0, 00) x (0, 00)) (31) 

denote the frequencies of, respectively, Si, S 2 and SK In terms of the Levy copula, A" 
takes the form 

A 11 = C(Ai, A 2 ). (32) 
The survival function F{~ of AS± is defined as 



F(Yj g ( Xl ,oo) x {0}) = ^((^, 00)) 
P(F 4 G(0,oo)x{0}) " A^ 



F l (xi) = F(AS 1 e [xx, 00)) = TO/>vr - fn _^ - rn1x = — , (33) 



where < X\ < 00. In terms of the Levy copula, F^ takes the form 

Fffa) = ^Fx(xi) - ^(AiAOri), A 2 ), (34) 

where Fx denotes the survival function of AS±. The distribution function F^~ of ASi is 
related to F^ by 

Ft{x x ) = P(A5j L G [0,x 1 ]) = 1 - Ff( Xl ). (35) 

Similarly, the distribution function Fx of ASx is given by Fx = 1 — A- At this point, 
we have related F^~ to C, Xx, A 2 and Fx. In the same way, the distribution function F^~ 
of AS 2 is related to C, Ai, A 2 and the distribution function F 2 of AS" 2 . In terms of the 
survival functions 

= 1 - F 2 X and F 2 = 1 - F 2 , (36) 

the relation takes the form 

F2M = ^F 2 (x 2 ) - ^C(\x, A 2 F 2 (x 2 )). (37) 
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Finally, the joint survival function F" of AS" is defined as 

F ll ( Xl ,x 2 ) =P(AS[' G (ari,oo),A5Jj G (x 2 ,oo)) 

_ f(Yj G (xi,oo) x (x 2 ,oo)) _ ^((a^oo) x (x 2 ,oo)) ( 38 ) 

~ p(r 4 g (o,oo) 2 ) ~ Ail ' 

where < ^i,X2 < oo. In terms of the Levy copula, F" takes the form 

F\x u x 2 ) = ^C(X 1 F 1 (x 1 ),X 2 F 2 (x 2 )). (39) 
The distribution function F" of AS 1 " is related to F" by 

F^{x!,x 2 ) = 1-F ll (x 1 ,0)-F ll (0,x 2 ) + F l, (^i,^2). (40) 

3.4 Examples of implied dependence 

A distributional survival copula C of F" satisfies 

fII( X1;X2 ) = C(fJ I (x 1 ),f|(x 2 )) (41) 

where 

Ff(x 1 ) = FW(x u 0) and Fj (x 2 ) = F 11 (0, x 2 ). (42) 

We assume here that Fi and F 2 are continuous, which implies that C is unique. For the 



Clayton Levy copula (29), substituting 



F i(xi) = ( ; — s = «i and 



2 

K° + a 2 5 / 

(A 2 F 2 (x 2 ))- 5 + Ar 5 ^ 1/5 



in the left- and right-hand side of Eq. (41) gives 



C(u 1 ,u 2 ) = (u^ 5 + u 2 s -l)- 1 / 5 , (44) 

which is the distributional Clayton copula. The distributional copula C of F" takes the 
form 

C(ux, u 2 ) = ((1 - Ul )- 5 + (1 - u 2 y 5 - iy 1/s + Ul + u 2 -l, (45) 

which collapses to uv for 5 4- and to min(u, t>) for 5 — >■ oo. The frequency A" implied by 
the Clayton Levy copula takes the form 

A" = (Ar 5 + A 2 5 )- 1/5 i (46) 
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which collapses to zero for 5 1 and to min(Ai, A2) for 5 — > 00. In summary, for 5 4- 0, 
the Clayton Levy copula implies A" = and independent components of S", while, for 
5 — > 00, it implies A" = min(Ai, A2) and comonotonic components of 5". 

As a second example, we consider the pure common shock Levy copula defined by 



Avanzi et al. (2011) as 



C(«i, u 2 ) = SuiU 2 l max (ui,u2)<oo + Mil« 2 =oo + u 2 l Ul =oo, where < 5 < min ^J-, . 

2 (47) 

For the pure common shock Levy copula, substituting 

Pf(xi) = Ui and F 2 (x 2 ) = u 2 (48) 



in the left- and right-hand side of Eq. (41) gives 



C{ux,u 2 ) = uiu 2 , (49) 

which implies that the components of 5" are independent. The frequency A" implied by 
the pure common shock Levy copula takes the form 

A 11 = 5AiA 2 , (50) 

which equals zero if 5 = and min(Ai,A2) if 5 = min(l/Ai, I/A2). In summary, the 
dependence implied by the pure common shock Levy copula is between frequencies only. 



4 Observation scheme and maximum likelihood esti- 
mation 

We consider a sample of jumps of a positive bivariate compound Poisson process discretely 
observed up to time t = T. The observation scheme is such that all jump sizes are 
observed, but it is not known which jumps stem from common shocks. We consider a 
partition of [0, T] in M intervals of equal length. The partition is chosen such that jumps 
of separate intervals can realistically be assumed not to stem from common shocks. In 
the context of operational risk modelling, with T/M either being a month or a quarter, 
this is the observation scheme typically assumed in the advanced measurement approach. 

The objective of this work is to estimate the parameters of the Levy copula in the 
observation scheme described above. A possible solution is to construct a likelihood 
function based on all possible combinations of jumps within each interval. If, within a 
certain interval, there are ri\ jumps within loss category one and n 2 jumps within loss 
category two, one can distinguish between min^!,?^) + 1 possibilities for the number of 



10 



common jumps < < rnin(rii, 712). Given a certain n\ there are 

(max(ni, ^2))! (min(ni, n 2 ))! 
(max(ni, n 2 ) — nl')! (min(n 1? n 2 ) — nil)! 

possibilities of distributing the common jumps over the observed jump sizes. Due to the 
large number of possibilities, a likelihood function based on all combinations of jumps 
is not feasible. An alternative approach is to construct a likelihood function based on 
the number of jumps and the expected jump sizes within the intervals. This approach, 
however, is also not feasible because the convolutions involved typically have no closed- 
form expressions. In the method proposed here, we use a sample consisting of the number 
of jumps and the maximum jump sizes within the intervals. For such a sample, we derive 
a closed-form likelihood function. Alternatively, a closed-form likelihood function based 
on the minimum jump sizes can also derived. In the context of operational risk modelling, 
however, one can expect the likelihood function based on maximum losses to me more 
variable with respect to model parameters than the likelihood function based on minimum 
losses. 



4.1 Discrete processes for maximum likelihood estimation 

We consider a partition II = {t , ti, .., £m} of [0, T], where t = 0, tM = T and tj — ti-i = 
T/M for all i = 1, . . . , M. The discrete process Z = (Zi)i=i,...,M = (Za, Zj 2 )i=i,...,M is 
defined as 



where 



Z^ = max ASHs) and Z%= max ASHs). (52) 



The random vector Zj is independent of and has the same distribution for all %. 

In the same way as Z, we define the discrete processes Z 1 - and Z" based on and zfj, 
respectively. These processes have the same properties as Z, but, in contrast to Z, they 
cannot be observed. The processes Z 1 - and Z" are independent. Also, the components of 
Z 1 - are independent. 

The discrete process iV = (Ni)i=i,...,M = (-^1,^2)1=1,...,]^ is defined as 

Nij = Nj(ti) - N^), (53) 

where the continuous process Nj counts the number of jumps of Sj. The random variable 
Nij has decomposition 

N%j = N± + JV|, (54) 
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where 

W/fe-i) and ivj =iV»(t i )-JV ll (*i-i)- (55) 

Here, the continuous processes iVj 1 and iV" count the number of jumps of Sj 1 and S", re- 
spectively. The random vector is independent of J r (t i _i) and has the same distribution 
for all i. In the same way as N, we define the discrete process N x based on N^. The 
process iV" = (ivf )j=i,... i M is the discrete process corresponding to iV'l The processes N 1 
and iV" have the same properties as N, but, in contrast to N, they cannot be observed. 
The processes N 1 - and iV" are independent. Also, the components of N 1 - are independent. 



4.2 The likelihood function 

Realizations of the process Z are collected in an M x 2 matrix z such that ^ is a realization 
of Zij. Similarly, the M x 2 matrix n holds realizations of iV. The likelihood function 
L Z) n corresponding to z and n takes the form 

M 

L z ,h — W^Zi 1 ,Zi2,n il ,hi2i (56) 
i=l 

where £z il ,z i 2,h il ,h t 2 denotes the likelihood function of the i-th row of z and n. The entries 
of L and C consist of Ai, A2, the parameters of Fi, the parameters of F 2 and, finally, the 



parameters of the Levy copula C. The function £ in Eq. (56) has no discrete time label 



i because the distributions of the random vectors do not depend on time. The likelihood 



function L factorizes cf. (56) because random vectors of different discrete time points are 



independent. The function C takes the form 

dF k Ax,oo) 



£x,y,k,l — l{fc=0,«=0}-^0,o(oO, (X)) + l{fc>CM=0}- 



dx 



(57) 



&F ,/(oo,2/) d 2 F k j(x,y) 

+ 1 {fe=o,i>o} ^ r -L{fc>cM>o} n n • 

ay oxoy 

Here, the function F kj i is defined as 

F ktl {x, y) =P(Z n < x, Z 12 < y, N n = k, N 12 = I) 

min(fc,i) 

= ^(Z n <x,Z 12 <y\N^ = k-J,N^ = l-n^Nf=n ]l )x, (58) 

nll=0 

P(iVi L 1 = k-J,N^ 2 = l-J,N ]l =J), 
where we have used that 

{N n = k,N 12 = l} = uX ( fc '°{A\i = k-n\N^ 2 = I - n",^" = J}. (59) 
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The functions Fk,o and Fqj with k, I > are assumed to be different iable with respect 
to the first and second entry, respectively. The functions F^i with k, I > are assumed 
to be continuously differentiable of second order. As we will see later, these assumptions 
hold if F\ and F 2 are continuously differentiable and the Levy copula C is continuously 
differentiable of second order on (0, Ai) x (0, A 2 ). 

To find an explicit expression for C, we need to calculate P(iVf 1 [ = k — n^,N^ 2 = 



= n") and to calculate the conditional probability appearing in Eq. (p8J). 
Because N x \, N X2 and ivf are statistically independent Poisson random variables with 
frequencies X^T/M, X 2 T/M and A"T/M, respectively, we find 



P(A> 1 - L 1 = k - n 11 , A\i = I - n 11 , iVj 1 = n 11 ) = 

(A^T/M) fc -""e-^ T / M {X^T/M) l -^e-^ T / M (A»T '/ 'M)" [l e ->"r/M 



• (60) 



(A; -nil)! (/-rill)! nil! 

The conditional probability appearing in i^o is given by 

P(Zn < x, Z 12 < 2/|A>n = 0, N 12 = 0, ivf = 0) = 1. (61) 

The conditional probability appearing in Fk,o with k > takes the form 

F{Z U < x, Z 12 < y\N^ = k, N^ 2 = 0, ivf = 0) = P(Z 1 ± 1 < x|i\^ = jfe) (62) 

where we have used 

{Z n < x,Z 12 <y}n {Nti = k,N^ 2 = 0, N 1 } = 0} = {Z^ < x,N& = k,N& = 0, A>} = 0} 

(63) 

and the independence of N± 2 and iVf from Zj^ and iV^. In terms of the distribution 



function F{~ of ASf defined by Eqs. (34) and (35), Eq. (62) becomes 



¥(Z U < x, Z 12 < y\N^ = k, = 0, N[ l = 0) = (F^(x)) h , (64) 

where we have used that the jump sizes of are iid. The conditional probability appear- 
ing in Fqj with I > is treated similarly. Finally, we consider a conditional probability 
appearing in Fj^i with k,l > 0. Making use of the independence of the processes , S 2 
and S", we find 

F(Z n < x, Z 12 < y\N^ = k- n ll ,N^ 2 = I - n 11 , JV"| = n") 

= P(Z 1 J 1 < x, Z\ x < x, < y, Z\ 2 < y\N^ — k — n", JV£ = I - n\ n\ = J) 

ii ii ~n (65) 

= P(Z 1 ± 1 < x\N^ — k- nl | )P(Z 1 J 2 < y\N X2 = I - n^P^j'i < x, Z\ 2 < y\N l l = n 11 ) 

= {Ft(x))^{Et{y)) l -«\F\z,y)f, 
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where F 2 is the distribution function of AS^ defined by Eqs. (36) and (37), and F" is 



the distribution function of AS"" defined by Eqs. (39) and MOl. 



Substituting the conditional probabilities of Eqs. ( |6~T|64|65[ ) in Eq. (58) and subse- 
quently substituting the resulting expression of F^i in Eq. fl57j), gives 



C x , y , k ,i = l{fc=o,i=o}P(iVi = 0, Nf 2 = 0, ivj = 0) 

+ l {fc >o,/=o}F(iV 1 - L 1 = k,N& = O.JVj = 0)fc(F 1 ± (x)) fc - 1 / 1 ± (x) 
+ l {fc =o,/>o}P(iV 1 - L 1 = 0,^ = l,Nf = o)Z(F 2 - L (j / )) , - 1 #(y) 

min(fc,i) 

+ l{fc>0,Z>0} 



(66) 



L1 = k- n\Nf 2 = I - n ll ,N{ = n")x 



nil =0 



d 2 



dxdy 



[(F 1 ± (x)) k - n \F^(y)) l - n \F^x,y)f } 



where f± and f 2 are the probability density functions corresponding to, respectively, F^~ 



and F 2 . From Eqs. (34) and (35), we find 



dFt(x) \ih{x) 



dx 



dC(u, A 2 ) 



du 



u=\iF±(x) 



(67) 



where fi is the probability density function corresponding to Fx and F\ = 1 — F\ is 
the survival function of A Si. The density fi exists because F\ and C are differentiate. 



Similarly, from Eqs. (36) and (37), we find 



dFi-{y) X 2 f 2 (y) 



dy 



1 - 



dC(X u v) 



dv 



v=\2F2{y) 



(6* 



where f 2 is the probability density function corresponding to F 2 and F 2 = 1 — F 2 is the 
survival function of AS 2 . The double derivative in the likelihood takes the form 

d 2 



dxdy 



(F 1 ± (x)) k - n \F 2 ± (y)) l - n \FHx,y)r 
= {k _J ){l _ n \\ )f ± {x)f ± {y){F ± {x)) ^^ 

+ n\l-n\\)f 2 \y)Fl{x,y){Ft{x)) k - n \F 2 \y)) l - n ^ 
+ n\k-n\\)ft{x)Fl{x^F 2 \y)) l - r \Ft{x)) k -^ 
+ n»(n« - l)Fhx,y)F 2 \x,y)(F 1 ± (x)) k - n \F 2 ± (y)Y- n \FHx,y)y 
+ n^f\x,y)(F^(x)) k - n \F 2 ± (y)) l - n \F\x,y)f-\ 



(69) 



\nll-2 



where 



F»(x,y) 



dFW(x,y) ^ 



dx 



F 2 \x : y) 



dFW(x,y) 
dy 



and f ll (x,y) 



d 2 FW(x,y) 
dxdy 



(70) 
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From Eqs. (39) and (40), we find 



F} (x,y) 



Ai/i(x) 
AH 



dC(u,X 2 F 2 (y))_ + dC(u,X 2 ) 



du 



du 



u=\±Fi (x) 



(71) 



Similarly, F 2 (x, y) takes the form 



F 2 0, y) 



A" 



dC(X 1 F l (x),vJ + dC(X 1 ,v) 



dv 



dv 



v=\2F2(y) 



(72) 



Finally, the density /'I is given by 

d 2 FW{ X) y) X 1 X 2 f 1 (x)f 2 (y) d 2 C(u,v) 



dxdy 



A" 



dudv 



(73) 



u=X 1 F 1 (x),v=\ 2 F 2 (y) 



The double derivative of Eq. (69) is continuous if F\ and F 2 are continuously different iable, 



and the Levy copula C is continuously differentiable of second order on (0, Ai) x (0, A2). 
The likelihood function C (and thereby L) is now completely specified in terms of the 
marginal frequencies, the marginal jump size distribution functions and the Levy copula. 
In the limit of continuous observation (M — >■ 00), the likelihood L converges to the 



likelihood derived by Esmaeili and Kluppelberg (2010) 



4.3 Maximizing the likelihood function 

The parameters of the process S can be estimated by maximizing the likelihood function 
L with respect to all its entries simultaneously. Alternatively, the likelihood function 
can be maximized with an approach similar to the inference functions for margins (IFM) 



approach of distributional copulas (Joe and Xu, 1996). The IFM method consists of two 
steps. In the first step, the parameters of Sj with j = 1,2 are estimated by maximizing 
the marginal likelihood function 



Kir 



Nj(T) 

(A;)^ (7 V A ' r n M>y)d, 

1=1 



(74) 



where Sj is the vector of jumps of Sj within [0,T]. In the second step, the estimates 
of the marginal parameters are substituted in L and the resulting likelihood function 
is maximized with respect to the Levy copula parameters. The IFM approach seems 
particularly suitable in the observation scheme of this work because the method makes 
use of all jump sizes (rather than the maximum jump sizes and the number of jumps in 
the M intervals) in estimating the marginal parameters. 
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5 A simulation study 



The quality of the estimation method of Section [4] is tested in a bootstrap analysis. The 
analysis consists of sampling many times from S on a period [0, 1] and estimating its 
parameters. The marginal jump size distribution function Fj with j — 1, 2 is given by 

F j {x) = l-e- B * (75) 



and dependence is introduced by the Clayton Levy copula of Eq. (29). The process S is 
thus described by five parameters (the marginal frequencies Ai and A2, the parameters 
61 and 82 of the marginal jump size distribution functions and the parameter 5 of the 
Levy copula). Given the parameters of the process, the sampling algorithm consists of 
the following steps: 

• Draw and from a Poisson distribution with frequency \± and A^ - , respec- 
tively. 

• Draw AT" from a Poisson distribution with frequency A". 

• Draw times from a uniform [0, 1] distribution. The resulting draws are the jump 
times of S±. The jump times of S£~ are determined similarly. 

• Draw iVf 1 times from a uniform [0, 1] distribution and apply the inverse of Fj 1 - to 
each draw. The resulting numbers are the jump sizes of S^. The jump sizes of 
are determined similarly. 

• Draw iV" times from a uniform [0, 1] distribution and apply the inverse of the 
marginal distribution function Ff defined as Ff(x) = F^(x, 00) to each draw. The 
resulting numbers Xj with i = 1, . . . , iV" are the jump sizes of Sf. (Note that the 
marginal distribution function F| defined here has one entry. In contrast, the func- 



tion Fj defined in Eq. ( 70 ) with two entries denotes the partial derivative of F" with 
respect to the first entry. We will use F| to denote both functions. The number of 
entries indicates to which function it refers.) 

Draw iV" times from a uniform [0, 1] distribution. The resulting draws are denoted 
by Ui with i = 1, . . . , N". Apply the inverse of the distribution function H Xi (y) = 
Ff(xi,y)/ff(xi), where fi(x) = dFf(x)/dx, to Uj for all i — 1, . . . , NK The resulting 
numbers yi with i — 1, . . . , are the jump sizes of S^. 
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In the last step, we have used that 



P(ASJj <y\ASl=x)= lim P(ASJj < y\x < ASf < x + Ax) 



lim 

AxlO 



II 
2 

x < AS l l <x + Ax, AS 2 l < y) 



N <f _|_ A nr.] 

(76) 



(x < AS} 1 < x + Ax) 



Um FW{x + Ax,y)-F\x,y) = Fj{x,y) = R 
A ™ Fj'(x + Ax) -FJ'(x) /J(x) J ' ! ' ' 



The algorithm to sample from S is the same as the one used by Esmaeili and Kliippelberg 



(2010) apart from the step where and are drawn. ( These authors first draw N±, 
N 2 and and then calculate and N% as = N x - iV" and A^ 1 = N 2 - NK This 
method may accidentally work, but it could result in negative N± and N 2 .) To prepare 
a sample to which we can apply the methods of Section |4j we perform the following steps 

• Divide the observation period [0, 1] in M intervals of equal length and determine Zn 
and Zi2 for alH = 1, . . . , M. This results in an M x 2 matrix z of maximum jump 
sizes. Also determine Nn and N i2 for all i — 1, . . . , M. This results in an M x 2 
matrix n of number of jumps. 

• Determine the vector s\ of all jump sizes of S% on [0, 1]. Similarly, determine the 
vector s 2 of all jump sizes of S 2 . (The vectors si and s 2 are used in the IFM 
approach.) 

The algorithm described above is repeated many times and based on each z, n, s± and 
s 2 , the parameters of S are estimated. Statistics of the resulting bootstrap estimates are 
given in Tables |T] and |TT] for the IFM approach and in Table III for maximum likelihood 



estimation. From Table [TJ we find that the difference between the bootstrap mean and 
the true value of 5 is within a standard error of the mean (bootstrap standard deviation 
divided by 10). This indicates that the estimate of 5 with the IFM approach is unbiased 
for the process under study and M = 100. From Table [TTJ, we find that the estimate of 5 
remains unbiased for other values of 5 and M. Also, for a fixed value of M, the bootstrap 
standard deviation is seen to be approximately proportional to 5, which means that the 



relative precision of the estimate of S does not depend on 5. From Table III , we find that 
the difference between the bootstrap mean and the true value of 5 is within a standard 
error of the mean (bootstrap standard deviation divided by vSO)- This indicates that the 
maximum likelihood estimate of 5 is unbiased for the process under study and M = 100. 
Similarly, from Tables [I] and III, estimates of the parameters of the marginal processes 



are seen to be unbiased in both methods. In terms of the bootstrap standard deviations, 
however, the IFM approach provides estimates of a slightly better quality. 
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Table I: Parameter value, bootstrap mean and bootstrap standard deviation of the param- 
eters of a bivariate compound Poisson process. The parameters consist of the marginal 
frequencies Ai and A2, the parameters 81 and 82 of the exponential marginal jump size 
distribution functions, and the parameter 6 of the Clayton Levy copula. The simulation 
period [0, 1] is divided in 100 intervals of equal length and the parameters are estimated 
with the IFM approach in 100 bootstrap samples. 





Ai 


A 2 


61 


82 


6 


value 


1000 


1000 


1 


1 


1 


bootstrap mean 


1004 


999 


0.999 


1.002 


1.007 


bootstrap standard deviation 


33 


36 


0.032 


0.031 


0.114 



Table II: Parameter value, bootstrap mean and bootstrap standard deviation of the Clay- 
ton Levy copula parameter 5 of a bivariate compound Poisson process. The process has 
marginal frequencies Ai = A2 = 1000 and marginal exponential jump size distribution 
functions with parameters 8\ = 82 = 1. The simulation period [0,1] is divided in M 
intervals of equal length and the parameter 5 is estimated with the IFM approach in 100 
bootstrap samples. For 5—1, the same bootstrap samples are used as in Table |TJ 



M 6 

value bootstrap mean bootstrap standard deviation 
50 1 0.989 0.143 

50 5 4.999 0.543 

100 1 1.007 0.114 

100 5 5.006 0.417 



Table III: Parameter value, bootstrap mean and bootstrap standard deviation of the 
parameters of a bivariate compound Poisson process. The parameters consist of the 
marginal frequencies Ai and A2, the parameters 8\ and 82 of the exponential marginal 
jump size distribution functions, and the parameter 5 of the Clayton Levy copula. The 
simulation period [0, 1] is divided in 100 intervals of equal length and the parameters are 
estimated with maximum likelihood estimation in 50 bootstrap samples. 





Ai 


A 2 


81 


82 


6 


value 


1000 


1000 


1 


1 


1 


bootstrap mean 


1004 


997 


0.993 


0.994 


0.992 


bootstrap standard deviation 


37 


39 


0.036 


0.042 


0.113 



6 A real data analysis 

We apply the methodology developed in Section [4] to the Danish fire loss data set publicly 
available at |http : //www . ma . hw . ac . uk/~mcneil| The data set consists of fire insurance 
data divided into loss of building, loss of content and loss of profit. Common shocks are 
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known in this data set and a Levy copula has already been fitted to the data based on 



the likelihood function of continuous observation (Esmaeili and Kliippelberg, 2010). To 
mimick an actual loss data set faced by banks in operational risk modelling, we consider 
a monthly observation and remove common shock information within months. 

6.1 Description of the data and pre-processing 

Details about the loss data set are described at its source http : / / www . ma . hw . ac . uk/| 



mcneil The data set consists of 2167 fire loss events over the period 1980 up to 1990 
(11 years). They have been adjusted for inflation and are given with respect to the year 
1985 in millions of Danish Kroner. Each loss event is divided into a loss of building, 
a loss of content and a loss of profit. In order to make a comparison with the case of 



known common shocks studied by Esmaeili and Kliippelberg (2010), we apply the same 



pre-processing to the data. This means that the loss of profit is not taken into account 
because it rarely has a non-zero value. For the remaining two loss categories, we consider 
only losses exceeding a million Kroner and take the logarithm of these losses. The resulting 
data set consists of 940 transformed loss events. 

To prepare a sample for the discretely observed process under study in this work, we 
consider a monthly partition of the 11 years and determine the maximum loss and the 
number of losses for each month and each loss category. This results in a 132 x 2 matrix 
z of maximum losses and a 132 x 2 matrix n of number of losses. Details about the 
maximum losses and the number of losses are given in Figure 1. To prepare for the IFM 
approach, we construct the vectors S\ and s 2 holding, respectively, all losses of building 
and all losses of content. The vector si holds 782 losses and the vector s 2 holds 456 losses. 

6.2 Marginal processes 

We will estimate a Levy copula with the IFM approach. This means that the parameters 
of the process Si of loss of building and the parameters of the process S 2 of loss of content 
are based on, respectively, si and s 2 - We use the same marginal jump size distributions 
as Esmaeili and Kliippelberg (2010). These authors use a Weibull distribution function 
for the log-losses and study the quality of the fit in terms of QQ-plots. The Weibull 
distribution function Fj of ASj takes the form 

Fj(x) = 1 - e-W a ^ j (77) 

where x > 0. By maximizing the marginal likelihood function 'Hj,s j , we find estimates of 



Xj, atj and /L- for j = 1, 2. The estimates are listed in Table IV 
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Figure 1: Monthly maximum log-losses and monthly number of losses for the Danish fire 
loss data set. The maximum log- losses are given as a function of time in the upper left 
panel (loss of building) and the middle left panel (loss of content). Similarly, the number 
of losses are given as a function of time in the upper right panel (loss of building) and the 
middle right panel (loss of content). The maximum log- losses of the two loss categories 
are plotted against each other in the lower left panel. Similarly, the number of losses of 
the two loss categories are plotted against each other in the lower right panel. 
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Table IV: Parameter estimates of the compound Poisson process corresponding to loss of 
building and the compound Poisson process corresponding to loss of content of the Danish 
fire loss data set. The annual frequency of losses of building is denoted by Ai and the 
annual frequency of losses of content by A 2 . The parameters ot\ and j5\ correspond to the 
Weibull distribution function of the log-losses of building. Similarly, the parameters a 2 
and 02 correspond to the Weibull distribution function of the log- losses of content. 



Ai 


A 2 at! 


01 


a>2 


ft 


estimate 71.1 


41.5 0.818 


1.197 


1.036 


1.131 



6.3 A goodness of fit test and selection of a Levy copula 

In the case of continuous observation or known common shocks, one can check if a Levy 
copula provides a reasonable fit to the data by inspecting the goodness of fit of the 
implied distributional copula between the components of AS 1 " . For a Clayton Levy copula, 



for example, the implied distributional copula is given by Eq. (45). The goodness of 



fit of a distributional copula can be assessed by transforming the correponding pseudo 



sample of probabilities to another sample of probabilities based on the copula (Breymann 



et al. 2010). If the copula provides a reasonable fit, the transformed probabilities are 
realizations of statistically independent uniform [0, 1] random variables. This can be 
tested by applying the inverse of the standard normal distribution function to the [0, 1] 
random variables and subsequently performing standard statistical tests. 

In the discrete observation scheme discussed in this work, common shocks are un- 
known. In order to test the fit and select a reasonable Levy copula, we construct a 
method similar to the goodness of fit test of distributional copulas. For this purpose, we 
define the distribution function Gk,i as 

G k>l (x,y) = F(Z n < x,Z l2 < y\N u = k,N 12 = I) = J k ^ X ' y \ . (78) 

For the i-th time interval, the probability Un is calculated as 

un = Gfi ai fi i2 (za, oo). (79) 
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We calculate the distribution function H Xj k,i{y) as 



H x , k ,i{y) = P(Zi2 < y\Z u = x, N u = k, N 12 = I) 
¥(Z n = x,Z 12 <y,N n = k,N 12 = l) 



u = x } N u = k,N 2 = I) 

,. P(x < Z n < x + Ax, Z i2 < y, N n = k, N 12 = I) (80) 

lim ~ = 

Ax+o p( x < z u < x + Ax, N n = k, N 12 = I) 

9F k ^(x,y) 
dx 



8F ktl (x,oo) ' 



dx 

The probability Ui 2 is calculated as 

Ui2 = H Ziuflll ,h l2 (z i2 ). (81) 

We select the rows of u that correspond to the rows of n of which both elements are non- 
zero. These rows correspond to the time intervals in which at least one loss is recorded 
in both loss categories. The selected rows are collected in a matrix v with two columns. 
In case of the Danish fire loss data, v is a 128 x 2 matrix. If the data generating process 
is correctly specified, v should be a realization of a random matrix with independent and 
uniform [0, 1] elements. To test this assumption, we translate v into a matrix w by 



w 



(82) 



where $ denotes the standard normal distribution function. The matrix w should have 
independent standard normal elements. For the Danish fire loss data, this is tested for 
the pure common shock Levy copula and the Clayton Levy copula. The results are listed 
in Table |V| For the pure common shock Levy copula, the estimated correlation coefficient 
P12 between the columns of w deviates from zero at the 0.01 level. This indicates that the 
pure common shock Levy copula is probably not appropriate for the data set. In contrast, 



the Clayton Levy copula provides a good fit. This is in line with the results of Esmaeili 



and Kliippelberg (2010). 



6.4 Estimation results 



Based on the analysis of Section 6.3, we use a Clayton Levy copula to fit the monthly 



Danish fire loss data. The parameters of the process are estimated with the IFM approach 



and the results are listed in Table VI We also estimate S with the IFM approach in case 
of known common shocks. This results in an estimate of 0.903, a bootstrap mean of 
0.904 and a bootstrap standard deviation of 0.043. These bootstrap statistics are based 
on 100 bootstrap samples. As expected, the bootstrap standard deviation with unknown 
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Table V: Test results of the goodness of fit of a bivariate compound Poisson process pa- 
rameterized with a Levy copula to the monthly Danish fire loss data set. The test is 
based on a 128 x 2 matrix w which, if the process is correctly specified, should contain re- 
alizations of statistically independent standard normal variables. The test results consist 
of the Jarque-Bera statistic (JB) of Jarque and Bera (1980), the mean fi, the standard 
deviation a and the serial correlation coefficient p of the two columns of w. A subscript 
indicates to which column the statistic refers. The test results also include the correlation 
coefficient pi2 between the columns of w. The p-values corresponding to the null hypoth- 
esis that the elements of w are realizations of independent standard normal variables, are 
given between parentheses below the corresponding statistics. The p-values are two-sided 
for fii, pi, p2 and pi2, and one-sided for &i and &2- 



Levy copula 


JB! 


JB 2 


h 


A2 


0-1 


0"2 


Pi 


h 


P12 


pure common shock 


0.79 


0.60 


-0.01 


-0.01 


1.04 


0.95 


0.09 


0.05 


0.25 




(0.64) 


(0.72) 


(0.94) 


(0.89) 


(0.26) 


(0.25) 


(0.29) 


(0.60) 


(0.01) 


Clayton 


1.14 


0.41 


-0.02 


-0.01 


1.00 


1.05 


0.08 


0.09 


-0.09 




(0.52) 


(0.80) 


(0.77) 


(0.92) 


(0.46) 


(0.21) 


(0.38) 


(0.31) 


(0.29) 



common shocks is larger than with known common shocks. The IFM estimate of 0.903 is 



close to the maximum likelihood estimate of 0.953 reported by Esmaeili and Kfiippeiberg 



(2010). 



Table VI: Estimate, bootstrap mean and bootstrap standard deviation of the parameters 
of a bivariate compound Poisson process fitted to the monthly Danish fire loss data. 
The parameters consist of the marginal frequencies Ai (loss of building) and A2 (loss of 
content), the parameters a.\ and j3\ of the Weibull distribution of log-losses of building, 
the parameters a 2 and (3 2 of the Weibull distribution of log- losses of content, and the 
parameter 5 of the Clayton Levy copula. The parameters are estimated with the IFM 
approach in 100 bootstrap samples. 





Ai 


A 2 


a 1 




Ot2 


fa 


5 


estimate 


71.1 


41.5 


0.818 


1.197 


1.036 


1.131 


0.695 


bootstrap mean 


70.5 


41.3 


0.818 


1.206 


1.038 


1.141 


0.699 


bootstrap standard deviation 


2.6 


2.1 


0.024 


0.034 


0.047 


0.046 


0.092 



7 Conclusions 

In summary, we have developed a method to estimate a Levy copula of a bivariate com- 
pound Poisson process in case the process is observed discretely with knowledge about 
all jump sizes, but without knowledge of which jumps stem from common shocks. The 
method is tested in a simulation study with a Clayton Levy copula. The results indicate 
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that the method is unbiased in small samples and that the bootstrap standard deviation of 
the Clayton Levy copula parameter is approximately proportional to its bootstrap mean. 
A goodness of fit test for the Levy copula is developed and applied to monthly log-losses of 
the Danish fire loss data set. The results indicate that the Clayton Levy copula provides 
a good fit to the data set. 

The method developed in this work is particularly useful in the context of operational 
risk modelling in which common shocks are typically unknown. To model dependencies 
between operational losses of different loss categories, the common practice in the banking 
industry is to use a distributional copula between either the number of losses or the ag- 
gregate losses within a certain time window. A disadvantage of this approach is that the 
distributional copula depends non-trivially on the length of the time window. If one has, 
for example, estimated a distributional copula between monthly losses, the distributional 
copula between yearly losses is typically unknown. A second disadvantage of the approach 
is that the nature of the model depends on the level of granularity. If one combines, for 
example, two loss categories connected by a distributional copula, the new loss category 
is typically not compound Poisson. These two issues are resolved by a multivariate com- 
pound Poisson process, which can be parsimoniously modelled with a Levy copula in a 
bottom-up approach. 



A Connection between Levy measures 

In this Appendix, we relate the Levy measures and to the Levy measures ui, 

v<i and the Levy copula C. On a Borel set (x±, oo) with < x\ < oo, the Levy measure 
is given by 

^((xi, oo)) = u((x 1: oo) x {0}), (83) 

which is equivalent to 

Vi((x\, oo)) = v((x\, oo) x [0, oo)) — lim i/((#i, oo) x [rr 2 , oo)). (84) 

In terms of 

vi({x\,oo) = v((xi,oo) x [0,oo)) (85) 
and the Levy copula, v±{{xi, oo) takes the form 

Vi ((x i, oo)) = ui((xi, oo)) — lim lim z/Qxi, oo), [x 2 , oo)) 

x'^xi X2IO 

= z/i((xi, 00)) - lim \imC(isi([x[, 00), v 2 ([x 2 , 00)) (86) 

x'^xi X2IO 

= ui((x 1 , 00)) - C(ui((x 1 , 00), i/ 2 (0, 00)), 
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Similarly, on a Borel set (x 2 , °o) for < x 2 < oo, the measure z/^ takes the form 

^((£2, 00)) = ^((^2, 00)) - C(z/i(0, 00), z/ 2 ((^2, 00))), (87) 

Finally, on a Borel set (x±, 00) x (x 2 , 00) with < x±, x 2 < 00, the measure i/H is given by 

z/"((x 1 ,oo) x (x 2 ,oo)) = lim lim [0^,00) x |V 2 ,oo)) 

= lim lim 00), i/ 2 [x 2 , 00)) (88) 

= C(U 1 ((X 1 , OO), l/ 2 (^2, OO)). 
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